Lower Bounds on the Ground State Entropy of the Potts Antiferromagnet on Slabs of 

the Simple Cubic Lattice 
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We calculate rigorous lower bounds for the ground state degeneracy per site, W, of the g-state 
Potts antiferromagnet on slabs of the simple cubic lattice that are infinite in two directions and finite 
in the third and that thus interpolate between the square (sq) and simple cubic (sc) lattices. We 
give a comparison with large-q series expansions for the sq and sc lattices and also present numerical 
comparisons. 

PACS numbers: 



I. INTRODUCTION 

Nonzero ground state entropy (per lattice site), So ^ 0, 
is an important subject in statistical mechanics, as an 
exception to the third law of thermodynamics and a phe- 
nomenon involving large disorder even at zero tempera- 

1 /n 

ture. Since Sq — ks hi W , where W = lim^oo VT tot , and 
n denotes the number of lattice sites, So is equiv- 
alent to W > 1, i.e., a total ground state degeneracy 
Wtot. that grows exponentially rapidly as a function of 
n. One physical example is provided by H2O ice, for 
which the residual entropy per site (at 1 atm. pressure) 
is measured to be So = (0.41 ± 0.03)/cb, or equivalently, 
W = 1.51 ± 0.05 A salient property of ice is that 

the ground state entropy occurs without frustration; i.e., 
each of the ground state configurations of the hydrogen 
atoms on the bonds between oxygen atoms minimizes the 
internal energy of the crystal [J]. 

A model that also exhibits ground state entropy with- 
out frustration and hence provides a useful framework in 
which to study the properties of this phenomenon is the 
g-state Potts antiferromagnet (5j-[7| on a given lattice A 
or, more generally, a graph G, for sufficiently large q. 
Consider a graph G = (V, E), defined by its vertex (site) 
and edge (bond) sets V and E. Denote the cardinalities 
of these sets as n(G) = \V\ = n and e(G) = \E\, and 



let {G} = lim 



n(ay 



, G. An important connection with 



graph theory is the fact that the zero-temperature par- 
tition function of the g-state Potts antiferromagnet on 
the graph G satisfies Z(G 7 q,T = 0) = P(G,q), where 
P(G, q) is the chromatic polynomial expressing the num- 
ber of ways of coloring the vertices of G with q colors 
such that no two adjacent vertices have the same color 
(called a proper q- coloring of G) Thus, 



W({G},q) = lim P(G,g) 



l/r, 



(1.1) 



In general, for certain special values of g, denoted g s , one 
has the following noncommutativity of limits [l(| 

lim lim P(G,q) 1/n ^ lim lim P{G,q) 1/n , (1.2) 

n— >-oo q— >q s Q^Qs Ti^oo 

and hence it is necessary to specify which order of limits 



that one takes in defining V7({G}, q). Here by ^({G}, q) 
we mean the function obtained by setting q to the given 
value first and then taking n — > 00. For the n — > 00 limit 
of a bipartite graph Gu P . , an elementary lower bound is 
W({Gbi P .}, q) > y/q — 1, so that for q > 2, the Potts anti- 
ferromagnet has a nonzero ground state entropy on such 
a lattice. A better lower bound for the square lattice is 
W7sg,g) > (g 2 -3g+3)/(g-l) [TTJ. In previous work - 
111 one of us and Tsai derived lower and upper bounds on 
W for a variety of different two-dimensional lattices. It 
was found that these lower bounds are quite close to the 
actual values as determined with reasonably good accu- 
racy from large-g series expansions and/or Monte Carlo 
measurements . 

In the present paper we generalize these lower bounds 
on two-dimensional lattice graphs by deriving lower 
bounds on W({G},q) for sections of a three-dimensional 
lattice, namely the simple cubic lattice, which are of infi- 
nite extent in two directions (taken to lie along the x and 
y axes) and finite in the third direction, z. By compar- 
ison with large-g expansions and numerical evaluations, 
we show how the lower bounds for the W functions for 
these slabs interpolate between the values for the (re- 
spective thermodynamic limits of the) square and simple 
cubic lattices. These bounds are of interest partly be- 
cause one does not know the exact functions W(sq, g) or 
W(sc, q) for general q. 



II. CALCULATIONAL METHOD 



Let us consider a section (slab) of the simple cubic 
lattice of dimensions L x x L y x L z vertices, which we de- 
note sc[(L x )bcx x (L y ) B cy x (L z )bc z ], where the bound- 
ary conditions (BC) in each direction are indicated by 
the subscripts. The chromatic polynomial of this lattice 
will be denoted P(sc[(L x ) B cx x {L v )bc v X (L z )BCy],g)- 
We will calculate lower bounds for W(sc[(L x )bcx X 
{L y ) B Cy x (L z ) B Cz],q) in the limit L x -> 00 and L y -> 00 
with L z fixed. These are independent of the bound- 
ary conditions imposed in the directions in which the 
slab is of infinite extent, and hence, for brevity of nota- 
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tion, we will denote the limit rim^ l v ^oo sc[(L x )bcx x 
{L y )BCy x (L z )bcz] simply as S (Lz)bc ^ where S stands 
for "slab". We will consider both free (F) and periodic 
(P) boundary conditions in the z direction, and thus 
slabs such as Ss F , S^p, etc. For technical reasons (to 
get an expression involving a trace of a coloring ma- 



To derive a lower bound on W(S(r,.) a ^_ , q), we gen- 
eralize the method of Refs. from two to three 
dimensions. We consider two adjacent transverse slices 
of the slab orthogonal to the x direction, with x values Xq 
and xo + 1 . These are thus sections of the square lattice of 
dimension L y xL z , which we denote G XQ — sq[{L y )Bc y x 
(L z )bcJx and G x „+i = sq[(L y ) B c y x (L z )bc.]x +1- We 
label a particular color assignment to the vertices of G Xf) 
that is a proper g-coloring of these vertices as C(G Xo ) 
and similarly for G XQ +\. The total number of proper 
q-colorings of G Xo is 

AT = P(G X0 ,q) = P(G X0+1 ,q) . (2.2) 

Now let us add the edges in the x direction that join 
these two adjacent transverse slices of the slab together. 
Among the Af 2 color configurations that yield proper q- 
colorings of these two separate yz transverse slices, some 
will continue to be proper (/-colorings after we add these 
edges that join them in the x direction, while others will 
not. We define an AT x A/"-dimensional coloring compat- 
ibility matrix T with entries Tc(g„ ),C(G x +i) equal to 
(i) 1 if the color assignments C{G Xo ) and C(G Xo+ i) are 
proper q-colorings after the edges in the x direction have 
been added joining G Xo and G Xa +i, i.e., if the color as- 
signed to each vertex v(xo,y, z) in G Xo is different from 
the color assigned to the vertex v{xq + 1, y, z) in G Xo +\; 
and (ii) if the color assignments C(G Xo ) and C(G Xo+ \) 
are not proper g-colorings after the edges in the x di- 
rection have been added, i.e., there exists some color as- 
signed to a vertex v(xq, y, z) in G Xo that is equal to a color 
assigned to the vertex v(xq + l,y, z) in G Xo +i. Clearly, 
Tij = Tji. The chromatic polynomial for the slab is then 
given by the trace 

P(sc[(L x ) P x (L y )BCy x (L z ) BG »],q) = Tr(T Lx ) . (2.3) 

Since T is a real symmetric matrix, there exists an or- 
thogonal matrix A that diagonalizes T: ATA -1 — Tdiag.- 
Let us denote the Af eigenvalues of T as At j , 1 < j < N- 
Since T is a real non-negative matrix, we can apply the 
generalized Perron-Frobenius theorem [13, [3 to infer 
that T has a real maximal eigenvalue, which we denote 



trix, as explained below) we will use periodic bound- 
ary conditions in the x direction. Note that the proper 
g-coloring constraint implies that FBC Z and PBC Z are 
equivalent if L z — 2. The number of vertices for G = 
sc[(L x )bcx x {Ly)Bc y x {L z ) B Cz] is n = L x L y L z . The 
specific form of Eq. (jl.ip for our calculation is 



(2.1) 

I 

X.T,max- It follows that 

lim [P(sc[(L x )px(L y ) B cyX(Lz)BC,},q)} 1/Ll = X T , ma x 

L x —}oo 

(2.4) 

Now for the transverse slices G XQ and G Xo +t, denoted 
generically as ts((L z )Bc z ): the chromatic polynomial has 
the form 

P(G X0) q) = P(G Xa+1 ,q) =^2 C J ( X t s «L*) BG ,),j) Lv 

(2.5) 

where the Cj are coefficients whose precise form is not 
needed here. The set of \ts((L z ) B a ),i' s 1S indepen- 
dent of the length L y and although this set depends on 
BC y , the maximal one (having the largest magnitude), 
^ts((L z ) BCz ),max, is independent of BC y (e.g., [3 and 
references therein). Hence, 

lim [P(G X01 q)]^ L y = lim (Af) 1 /^ 

Ly >-OG Ly" 7-CXj 

= \s{(L z )Bc z ),max • ( 2 -6) 

The two adjacent slices together with the edges in 
the x direction that join them constitute the graph 
sc[2p x {Ly)BCy x [L z )bc z ]- We denote the chromatic 
polynomial for this section (tube) of the sc lattice as 
P(sc[2 F x (Ly) B cy x (L z ) BC J,q) (which is equal to 
P(sc[2p x {Ly)scy x (L z )Bc z ],q) because of the proper 
g-coloring condition). This has the form 

P(sc[2 F x (L y ) B cy x (L z ) B c z ],q) 

= ^2 C '](^tube{(L z ) BCz ),j) Ly (2-7) 

j 

where are coefficients analogous to those in (|2.5p . 
Therefore, 

lim [P(sc[2 F x (Ly)BCy x (L z ) B c z },q)} 1/Ly = 

Ly— tOC 

— ^tube((L z ) B c z ),max ■ (2-8) 

Now let us denote the column sum (CS) 

Af 

CS j (T) = ^2T ij , (2.9) 

i=i 



I 

ff(S(i,) BCl , q) = lim lim [P{sc[{L x ) P x {L y ) B cy x (L z ) BCz ], q)} 1 

L v — too L x —>oo 
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which is equal to the row sum Y^jLi Tij, since T 1 = T. metric matrix T and r G N+ [li 
We also define the sum of all entries (SE) of T as 



SE(T) = ]T Ty 



(2.10) 



SE{T r ) 
Af 



1/r 



(2.12) 



Note that SE(T)/N is the average row (= column) sum. The lower bound is then 
Next, we observe that 



SE{T) = P(sc[2 F x (L y ) BCy x (L z ) BCz ],q) . (2.11) 

To obtain our lower bound, we then use the r = 1 
special case of the theorem that for a non-negative sym- where 

I 



W(S (Lt)BOt ,q)>W(S iLM)BOm ,q)t 



(2.13) 



lim 



P(sc[2 F x (Ly) BC y x (£ z )BcJ,g) 

P(sy[(L y ) S Cy X (L z ) BCz ],q) 



l/(L y L z 



^tube((L z ) bc z )>T7 la 3 ; 
Ats((L z ) bo z ) ,max 



l/L z 



(2.14) 



III. RESULTS FOR SLAB OF THICKNESS 
L z = 2 WITH FBC Z 



We now evaluate our general lower bound in Eqs. 
()2.13j) and (|2.14j) for a slab of the simple cubic lattice 
with thickness L z = 2 and FBC 2 , denoted S 2f . In this 
case the transverse slice is the graph sq[2 F x (L y ) B c y ]- 
For FBCj,, an elementary calculation yields 

P(sq[2 F X (L y ) F } 1 q) = q(q~l)(q 2 ~3q + 3) L y- 1 (3.1) 



with a single \ ts (2 F ) = 
the maximal A for PBC y [1 



^ = A ts (2), and this is also 
20], so that 



\s(2),max — q 2 — '&q + S ■ 



(3.2) 



be 



A 



tube{2) ,max 



8q 3 + 29q 2 - 55y + 46 



where 



'R 



22 



R 22 = q 8 - 16q 7 + 118g 6 - 526 9 5 + 1569g 4 
- 3250g 3 + 4617g 2 - 4136g + 1776 . 



(3.4) 



(3.5) 
and 



We then substitute these results for A ts ( 2 ),max 
\ube(2),max into the L z — 2 special case of (|2.14[) to 
obtain W(S 2 ,q)e, and thus the resultant lower bound 
on W(S 2F ,q) = W(S 2p ,q) = W(S 2) q); W(S 2) q) > 
W{S 2 ,q) e . 



We next use the calculation of 

P(sc[2 F x (L y )p x 2 F ],q) = P(sc[2 F x 2 F x (L y ) F ],q) 

= P{sq[A P x (L y ) P ],q) 

(3.3) 

in Ref. [13] (where each of the 2 F BC's is equivalent to 
2p), from which we calculate the maximal Xtube(2).max to 



IV. COMPARISON WITH LARGE-g SERIES 
EXPANSIONS 

One way to elucidate how this lower bound W(S 2 ,q)e 
compares with the exact W(sq, q) and W(sc, q) is to com- 
pare the large-g series expansions for these three func- 
tions. For this purpose, it is first appropriate to give 
some relevant background on large-y series expansions for 
W({G},q) functions. Since there are q n possible color- 
ings of the vertices of an n-vertex graph G with q colors 
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if no conditions are imposed, an obvious upper bound 
on the number of proper g-colorings of the vertices of G 
is P{G,q) < q n . This yields the corresponding upper 
bound W({G},q) < q. Hence, it is natural to define a 
reduced function that has a finite limit as q — > oo, 



W r ({G},q)=q- 1 W({G},q) . 



(4.1) 



For a lattice or, more generally, a graph whose vertices 
have bounded degree, W r ({G}, q) is analytic about 1/q = 
0. (W r ({G}, q) is non-analytic at 1/q = for certain 
families of graphs that contain one or more vertices with 
unbounded degree as n — > oo, although the presence of 
a vertex with unbounded degree in this limit does not 
necessarily imply non-analyticity of W r ({G}, q) at 1/q = 
[H], [22j].) It is conventional to express the large-g Taylor 
series for a function that has some factors removed from 
W r , since this function yields a simpler expansion. A 
chromatic polynomial has the general form 



n-fc(G) 

P{G,q)= Y, (-l) i On-i«* i 



(4.2) 



where the a„_j > and k(G) is the number of connected 
components of G (taken here to be k(G) — 1 without 
loss of generality). One has a n = 1, a n -i = e(G), and, 
provided that the girth g(G) > 3 [23], as is the case here, 
a„_2 = ( 2 )■ ^ ^-regular graph is a graph such that 
each vertex has degree (coordination number) n. For a 
K-regular graph, e(G) = nn/2. The coefficients of the 
three terms of highest degree in q in P(G, q) for a k- 
regular graph are precisely the terms that would result 
from the expansion of [q(l — g ,-1 ) re / 2 ]". Hence, for a re- 
regular graph or lattice, one usually displays the large-g 
series expansions for the reduced function 



W(A,q) 



W(A,q) 
g(l-g-l)«/2 



(4.3) 



The large-g Taylor series for this function can be written 
in the form 



w{K,q) = i + Y w ^y j > 

3=1 



where 



y 



1 



(4.4) 



(4.5) 



The two results that we shall need here are the large- 
q (i.e., small-?/) Taylor series for W(sq, q) and W(sc,q). 
The large-g series for W(sq, q) was calculated to succes- 
sively higher orders in [2J]-[28(. Here we only quote the 
terms to 0(y n ): 



W(sq,q) = l + y 3 + y 7 + 3y 8 + 4y 9 + 3y 
+ 3y n +0(y 12 ) . 



As noted above, lower bounds on W(A, q) obtained 
from the inequality (|2.12l) for two-dimensional lattices A 
were found to be quite close to the actual values of the 
respective W(A, q) for a large range of values of q. This 
can be understood for large values of q from the fact that 
they coincide with the large-g expansions to many orders, 
and the agreement actually extends to values of q only 
moderately above q = 2. For example, the lower bound 
on W(sq,q) in fllj is equivalent to W(sq, q) > (1 + y 3 )- 
This agrees with the small-?/ series up to order 0(y e ), as 
is evident from comparison with Eq. (|4.6[) . This lower 
bound also agrees quite closely with the value of W(sq, q) 
determined by Monte Carlo simulations in [l(J Il3l| 
(see Table 1 of [l| and Table 1 of [H]). We include 
this comparison here in Table HI For our pur poses, it 
is sufficient to quote the results from Ref. [l(| only to 
three significant figures. Since we are using large- q series 
for this comparison, we list the results in Table U for 
a set of values q > 4. As another example, the lower 
bound obtained for the honeycomb lattice in Ref. [l3[, 
W(hc,q) > (1 + y 5 ) 1 / 2 , agreed with the small-?/ series 
for W(hc,q) to 0(y 10 ). Thus, it was found that for all 
of the cases studied, W(A, q)e provides not only a lower 
bound on W(A,q), but a rather good approximation to 
the latter function. It is thus reasonable to expect that 
this will also be true for the lower bounds W(S(l z ) bc , q)i 
for the slabs Sl z of the simple cubic lattice considered 
here, of infinite extent in the x and y directions and of 
thickness L z in the z direction. 

From ingredients given in Ref. [25[ , we have calculated 
a large-g expansion of the W(sc, q) for the simple cubic 
(sc) lattice and obtain 



W(sc, q) = 1 + 3j/ 3 + 22y 5 + 31y 6 + 0{y 7 ) 



(4.7) 



In Table U we list the corresponding values of W(sc,q) 
obtained from this large-g series, denoted W{sc, q) ser ., 
for q > 4. We also list estimates of W(sc, q), denoted 
W(sc,q)MC for 4 < q < 6 from the Monte Carlo calcu- 
lations in Ref. [2^] . One sees that the approximate values 
obtained from the large-q series are close to the estimates 
from Monte Carlo simulations even for q values as low as 
</ 1- 

The coordination number for the S2 F slab of the sim- 
ple cubic lattice (of infinite extent in the x and y direc- 
tions) is k(S2 f ) = 5. We thus analyze the reduced func- 
tion W(S 2 , q)t = W(S 2 ,q)t/[q(l - q- x f /2 ]. This has the 
large-g (small-y) expansion 



W(S 2 ,q)i = l + 2y 3 + 2y 5 + 9y 6 + 0(y 7 ) 



(4.8) 



(4.6) 



As this shows, W(S2,q)i provides an interpolation be- 
tween W(sq, q) and W(sc, q); for example, the coefficient 
of the y 3 term is 1 for A = sq, 2 for A = S2, and 3 for 
A = sc. Furthermore, the coefficient of the y 5 term is 
for A = sq, 2 for A = S 2 , and 22 for A = sc. This 
is in agreement with the fact that the exact functions 
W(S^ Lz - )F ,q) interpolate between W(sc,q) and W(sc,q) 
as L z increases from 1 to 00 (30j and the expectation, 
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as discussed above, that W(S^ ) f , q)i should be close to 
W(S {Lz)F ,q). 



RESULTS FOR SLABS OF THICKNESS 

L z = 3, 4 WITH FBC 2 



One then needs \tube(3 P ),max- The relevant (4x4 di- 
mensional) transfer matrix for this tube graph was cal- 
culated for Ref. (34[, and we have used this to compute 
\ube(3 P ),max numerically. The results for W(S 3p ,q) are 
given in Table [I] 



For the slab of the simple cubic lattice with thick- 
ness L z = 3 and FBC Z , denoted S3 F , the transverse 
slice is the graph sq[3p x (L y )BCy\- The chromatic poly- 
nomials P(sq[3p x (L y )F],q), P(sq{3F x (L y )p],q), and 
P(sq[3p x (L y )xp],q) (where TP denotes twisted peri- 
odic, i.e., Mobius BC) were computed for arbitrary L y 
in Refs. [3l|, (32|, and [33j], respectively, and the maxi- 
mal A was shown to be the same for all of these bound- 
ary conditions. For the reader's convenience, we list this 



\ s {z F ).max in Eqs. (|9.2p and (|9.3[) of the Appendix. The 
other input that is needed to obtain the lower bound in 
Eq. (|2.14|) is the maximal A for the chromatic polynomial 
of the sc[2 F x3 F xL y ] tube graph, i.e., \tube(z P ),max- The 
relevant transfer matrix that determines the chromatic 
polynomial for this tube graph was given with Ref. [3~i| . 
Because it is 13 x 13 dimensional, one cannot solve the 
corresponding characteristic polynomial analytically to 



obtain A 



culate A 



tube(3 F ) ,max 



for general q. However, one can cal- 
tube{3 F ),max numerically, and we have done this. 
Combining these results with Eqs. (|9.2[) and (I9.3[) . we 
then evaluate the lower bound W(S3 F ,q)i by evaluating 
the L z = 3 special case of (I2.14[) . 

For the slab of the simple cubic lattice with thickness 
L z = 4 and FBC Z , S4 F , the transverse slice is the graph 
sq[Ap x {Ly)BCy\- Here the maximal \ts(A F ),max is the so- 
lution of the cubic equation (|9.5p given in the Appendix. 
One also needs the maximal A for the chromatic polyno- 
mial of the sc[2 F x4 F xL y ] tube graph, i.e., \ tu be(4 F ), m ax- 
The relevant (136 x 136 dimensional) transfer matrix for 
this tube graph was calculated for Ref. [lij , and we have 
used this to compute Xtube(4 F ),max numerically. We then 
obtain the lower bound W{S>4, F , q)e from the L z = 4 spe- 
cial case of Eq. (|2.14[) . The results for W(S 3F ,q) and 
W(S 4f , q) are listed in Table H 



VI. RESULT FOR SLAB OF THICKNESS L z = 3 
WITH PBC Z 



It is also of interest to obtain a lower bound for W 
for a slab with periodic boundary conditions in the z 
direction, since these minimize finite-volume effects. For 
this purpose we consider the slab of the simple cubic 
lattice with thickness L z = 3 and PBC 2 , S3 P . In this 
case the transverse slice is the graph sq[3p x {L y )BCy\- 
For FBC a the chromatic polynomial involves only one A, 
and this is also the maximal A for PBCj, and TPBCj, [36| . 
viz., 

hs(3 P ),max = Q 3 ~ ^ + 14<Z - 13 . (6.1) 



VII. DISCUSSION 

Since the slabs of infinite extent in the x and y direc- 
tions and of finite thickness L z geometrically interpolate 
between the square and simple cubic lattices, it follows 
that the resultant W functions for these slabs interpo- 
late between W(sq,q) and W(sc,q) 30]. Given that it 
was shown previously that the lower bounds W(A, q)j> 
obtained by the coloring matrix method are quite close 
to the actual values of the respective W(A, q) for a num- 
ber of two-dimensional lattices, this is also expected to 
be true for the W(S( Lz ) BCr ,q)g bounds. We have shown 
above how W(S2,q)e interpolates between W(sq,q) and 
W(sc, q) via a comparison of the large-g series expan- 
sions for these three functions. Table fl] provides a fur- 
ther numerical comparison for W(Sa F , q)i, W(S4 F ,q), 
and W(S3 P , q) with W(sq, q) and W(sc, q), the latter be- 
ing determined to reasonably good accuracy from large- 
q series expansions and, where available, Monte Carlo 
measurements. As noted, the lower end of the range of 
q values for the comparison is chosen as q — 4 in view of 
the use of large-g series. 

For sections of lattices, and, more generally, graphs 
that are not K-regular, one can define an effective vertex 
degree (coordination number) as [lij 



c e// 



2e(G) 



(7.1) 



For 3 < L z < oo, the slab of the simple cubic lattice (of 
infinite extent in the x and y directions) with FBC Z is 
not K-regular, but has the effective coordination number 



K eff(S(L z 



2 3 



1 

T 2 



(7.2) 



We observe that for the q values considered in Table 
HI W(sq,q) > W(S 2 ,q) e > W(S3 F ,q)e > W(S iF ,q) > 
W(sc,q). The fact that for fixed q, the exact function 
W(S(la f , q) is a non-increasing function of L z , and, for 
q > 2, a monotonically decreasing function of L z , follows 
from a theorem proved in Ref. [30(. To the extent that 
the lower bounds W{Sil ! \ f , q)e lie close to the actual val- 
ues of W(S(Lz)f> q), it is understandable that they also 
exhibit the same strict monotonicity. As was noted in 
Ref. [30(, the reason for the monotonicity of the exact 
values is that the number of proper g-colorings per ver- 
tex of a lattice graph is more highly constrained as one 
increases the effective coordination number of the lattice 
section. (This is also evident in Fig. 5 of [10|.) In the 
present case, the monotonicity can be seen as a result of 



the fact that the effective coordination number increases 
monotonically as a function of L z . 

The use of periodic boundary conditions in the z direc- 
tion minimizes finite-size effects, so that for a given L z , 
W(S( L ^ p ,q) would be expected to be closer to W(sc,q) 
than W(S(l^ f , q) 30]. Again, to the extent that the 
lower bounds are close to the actual W functions for 
these respective slabs, one would expect W(S^ z ) p ,q)e 
to be closer than W{S( L ^ F , q)g to W(sc,q). Our results 
agree with this expectation. In contrast to W(S(l z ) f > <z)> 
W(S/l\ p , q) is not, in general, a non-increasing function 
of L z , as was discussed in general in [3(| (see Fig. 1 
therein). Thus, values oiW{StL z \ p ,q), and hence, a for- 
tiori, W(S(L m ) P , q)g, may actually lie slightly below those 
for W(sc,q), as is evident for the W(S3 P ,q)e entries in 
Table IB 



with 

R3 = (q 2 -5q + 7)(q i - 5q 3 + llq 2 - I2q + 8) . (9.3) 

For the infinite-length strip of the square lattice with 
width 4 and free transverse boundary conditions, sq[Ap x 

00] dH 

W{sq[A F X 00], q) = (X iF , max ) 1/4 (9.4) 
where \4 F . m ax is the largest root of the cubic equation 



with 



04 F ,lX + 4f ,2X + 4f ,3 



= 



(9.5) 



7q 3 - 23q 2 + Alq - 33 (9.6) 



VIII. CONCLUSIONS 

In this paper we have calculated rigorous lower bounds 
for the ground state degeneracy per site W, equivalent 
to the ground state entropy So = ks InW, of the g-state 
Potts antiferromagnet on slabs of the simple cubic lattice 
that are infinite in two directions and finite in the third. 
Via comparison with large-q expansions and numerical 
evaluations, we have shown how the results interpolate 
between the square (sq) and simple cubic (sc) lattices. 



°4 F ,2 

and 

&4 F ,3 



2q 6 - 23q 5 + 116g 4 - 329g 3 + 553g 2 - 517q + 207 

(9.7) 

-q s + 16q 7 - 112g 6 + 449g 5 - 1130<? 4 + 1829g 3 
1858g 2 + 1084g - 279 . (9.8) 
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IX. APPENDIX 



We note the following results on E d lattices and lat- 
tice sections: W{K.bi P . , 2) = 1 for any bipartite lattice; 
W(sq,3) = (4/3) 3 / 2 U; and W{{L},q) = W{{C},q) = 
q — 1, where L n and C n denote the n- vertex line and cir- 
cuit graphs. For the infinite-length square-lattice strip 
of width 2, W(sq[2 F x 00], q) = W(sq[2 P x 00], q) — 
\J q 2 — 3q + 3, where, as in the text, the subscripts F 
and P denote free and periodic boundary conditions in 
the direction in which the strip is finite. For the infinite- 
length strip of the square lattice with (transverse) width 
3 and free transverse boundary conditions, sq[^F X 00] 

W(sq[3 F X 00], q) = {\ 3F , max ) 1/3 (9.1) 

where 



*3 F : max 



(q - 2){q 2 - 3q + 5) + ^/R 3 



(9.2) 



7 



TABLE I: Comparison of lower bounds W(S(l z ) bc ,q)i for 
(Lz)bCz = 2p = 2p, 3p, Ap, 3p with approximate values of 
W(A, q) for the square (sq) and simple cubic (sc) lattices A, as de- 
termined from large-g series expansions, denoted W(A, q) S er. and, 
where available, Monte Carlo simulations, denoted W(A, q)mc- 
We also list W(sq, q)i for reference. See text for further details. 



q 


W(sq,q) M c 


W(sq, q) S er. 


W(sq,q) e 


W(S 2F ,q) e 


W{Sa F ,q)i 


W(S 4F ,q) e 


W(SC, q) aer . 


W(sc,q) 


mc W(S 3p ,q)e 


4 


2.34 


2.34 


2.33 


2.13 


2.07 


2.04 


2.06 


1.9 


1.78 


5 


3.25 


3.25 


3.25 


2.96 


2.875 


2.83 


2.75 


2.7 


2.62 


6 


4.20 


4.20 


4.20 


3.87 


3.765 


3.71 


3.58 


3.6 


3.51 


7 


5.17 


5.17 


5.17 


4.81 


4.69 


4.64 


4.48 




4.43 


8 


6.14 


6.14 


6.14 


5.76 


5.64 


5.58 


5.41 




5.37 


9 


7.125 


7.125 


7.125 


6.73 


6.605 


6.54 


6.36 




6.325 


10 


8.11 


8.11 


8.11 


7.71 


7.58 


7.51 


7.32 




7.29 


100 




98.0 


98.0 


97.5 


97.4 


97.3 


97.0 




97.0 
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